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Abstract: The conjecture is verified that the optimum, energy minimizing mag- 
netic flux for a half-filled band of electrons hopping on a planar, bipartite graph 
is Tv per square plaquette. We require only that the graph has periodicity in one 
direction and the result includes the hexagonal lattice (with flux per hexagon) 
as a special case. The theorem goes beyond previous conjectures in several ways: 
(1) It does not assume, a-priori, that all plaquettes have the same flux (as in 
Hofstadter's model); (2) A Hubbard type on-site interaction of any sign, as well 
as certain longer range interactions, can be included; (3) The conclusion holds for 
positive temperature as well as the ground state; (4) The results hold in D > 2 

ON 



O 

O 



dimensions if there is periodicity in D — 1 directions (e.g., the cubic lattice has 
the lowest energy if there is flux tc in each square face). 
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The flux phase conjecture states that the ground state (g.s.) energy minimizing mag- 
netic flux through a square, planar lattice on which free electrons hop is 7r per plaquette 
when the electron filling factor is 1/2 [1-4]. (Zeeman terms are excluded.) This conjec- 
ture, along with extensions to positive temperature, higher dimensional geometries and 
allowance for some electron-electron interactions, will be proved here. 

If the sites of the lattice are interpreted as atoms in a solid then flux tt would corre- 
spond to magnetic fields available only on neutron stars. The significance of the flux phase 
is thus not primarily as a literal interpretation in terms of physical magnetic fields. One 
interesting interpretation concerns mean field calculations connected with superconductiv- 
ity. The main interest, however, in the author's view, is that it shows that diamagnetism 
(which states that the optimal flux is zero — and which is correct when the electron den- 
sity is very small) can be reversed when the density is high. Indeed, it can be maximally 



1 



FLUX PHASE BHL14/Aug/94 



reversed, as in this case (since flux on a lattice is determined only modulo 2tt). Thus, 
there is a peculiar, non-intuitive and poorly understood effect of the Pauli principle on the 
way in which orbital motion interacts with magnetic fields. It has been studied extensively 
[5-13]. See [10] for some history. 

To define things precisely, we start with a general finite graph A, which is a collection 
of |A| sites and certain bonds denoted by xy with x and y in A and x 7^ y. A positive 
weight \t xy \ — \t yx \ is specified in advance for each bond. By convention t xx = 0. The 
hopping amplitude is then t xy = \t xy \ exp[i(j)(x, y)], with (j)(x, y) = — 4>(y, x) for hermiticity, 
and the problem is to find the numbers <p(x,y) that minimize the total electronic ground 
state energy (when (3 = 1/kT = 00) or free energy (when (3 < 00). 

A circuit in A is a sequence of points xi, X2, • • • , x n , x\ with t XiXi+1 7^ for all i. The 
flux through this circuit is Yl7=i 4>(xi, Xi + i) (mod27r). It is a fact [10] that the spectrum 
of the hermitian matrix T = {t xy } Xjye A depends on the 0's only through the fluxes. This 
is also true of the Hamiltonians below. No a priori assumption is made that the flux need 
be the same in all plaquettes; indeed, the flux is not even assumed to be the same for up- 
and down-spin electrons. We allow different | tail's and 4>{x, y)'s for the up- and down-spin 
electrons. We denote these by T\TK Thus, our results apply to the Falicov-Kimball 
model (where = 0), for example. 

The electronic kinetic energy operator, in second-quantized notation, is 

K = - E ^4 T ST+^4iSi- C 1 ) 

x,yeA 

The c's satisfy the fermion anticommutation relations {c xa , c yT } = 5 xy 5 aT , {c xcr , c yT } = 0. 
The electron number is N = N-\ + N± and N a = J2 x eA n ^ with n xu = cLc_. If the 
Hamiltonian H equals K, then the ground state energy £0 would be E = J2\<o H^ ) + 
X!a<o ^(T^), i- e -? t ne sum of t ne negative eigenvalues, A(T), of the matrices and TK 

For bipartite graphs (i.e., A = A U B, An B = and t xy = unless x G A, y G B or 
x G B, y G A) Eq is achieved when N = |A|, and hence the appelation half filled band, 
since < N < 2|A| in general. In the Hubbard model H = K + W°, 

W° = J2u x (n x i-±)(n xl -±), (2) 

xeA 

but we can also add certain longer range density-density interactions, W d and spin-spin 
interactions W s to be specified later, of the form (with w^. y = w^.* = w yx , k = d, 1,2,3 
and being Pauli matrices) 

W d = O-T + n xi - 1)(% T + n yl - 1) 

x,yeA 

^ S = EE E <(4«si T c XT )(cisi xCyX ). (3) 

j=l x,yeAa,T,fi,X 
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W d and W s are invariant under the unitary hole-particle (h-p) transformation r, 
with TC xcr r _1 = c\. a . The kinetic energy operator with complex T's satisfies 

tK(T\ T^r -1 = K(—T^*, —T- 1 *) (4) 

where T* denotes the complex conjugate matrix t xy — > t* y = t yx . 

The grand-canonical partition function of our system, which we want to maximize, is 
Z = Trexp[— f3H] at inverse temperature (3. By h-p symmetry, (N) = TrNe -1311 /Z is |A|, 
which is the half-filled band. The /3 — > oo limit is discussed at the end. 

Henceforth, all graphs will be bipartite, in which case all elementary circuits contain 
an even number of sites and bonds. The original flux phase conjecture is that when 
H = K, N = |A| and A is planar the optimum choice of fluxes is n in every circuit 
containing 0(mod4) sites and in circuits with 2(mod4) sites. 

Several cases of this were proved in [9, 10] and it was pointed out in [10] that the 
conjecture can not always hold for arbitrary values of \t xy \. It depends on A. Despite this 
caveat, however, it was proved in [10] that log [det(T 2 )] = ^2jlog[\j(T) 2 ] is maximized 
when the fluxes accord with the conjecture. This is true for an arbitrary bipartite, planar 
graph with arbitrary \t xy \'s. (On any bipartite graph, the nonzero eigenvalues of T come 
in opposite pairs A, —A.) 

The generalized flux phase conjecture is that the above choice is also optimal for 
H = K + W° + W d + W s and for all j3 < oo. We shall prove this here for graphs that 
have a certain periodicity. The usual square lattice with periodic boundary conditions is 
included. The result holds also for higher dimensional, non-planar graphs. The type of 
graphs A considered here is illustrated in Fig. 1 for the planar case. A is wrapped on a 
cylinder, i.e., the sites at the right end are identified with the sites on the left. The Idyl's 
on the vertical edges must be periodic, i.e., they are allowed to vary in an arbitrary way 
along each column, but all the columns must be identical. The horizontal Idyl's can also 
vary as we move vertically but they are only required to have period 2 in the horizontal 
direction, i.e., every second column of horizontal edges must be the same. Thus, if we 
erase the horizontal edges 1, 3, 5, etc. from every second column and if we erase edges 
2, 4, 6, etc. from the remaining columns in between, the hexagonal lattice is obtained. 
Our result includes this case, and flux tt in each (imaginary) square then implies flux 
in each hexagon. Octagons, decagons, etc. can also be included provided periodicity 2 is 
maintained. We can, if we wish, insert vertical edges connecting the bottom row to the 
top row, as indicated by the vertical arrows in Fig. 1. A will no longer be planar, but that 
does not matter. 

An easier way to state the periodicity requirements is to cut A at the two dashed lines, 
called P in Fig. 1. The two half-cylinders (as well as the Idyl's on the edges) are required 
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to be mirror images of each other. The (tail's on the edges that intersect the dashed lines 
are arbitrary because each is its own mirror image. We then say that \T\ is reflection 
symmetric with respect to the cutting lines, P, through the bonds. Our theorem says that 
when H = K the optimum flux is then 7r for those squares containing the cutting lines. If 
the |ta;j/|' s are reflection symmetric w.r.t. every choice of cutting lines - which is equivalent 
to the above periodicity requirement - then flux n will be optimal in every square of A. 

The periodicity of the Idyl's mentioned above is not needed for the theorem below. 
Only reflection symmetry is needed. The periodicity comes in when we wish to insure flux 
7T in every plaquette of A. This is achieved by repeated reflection in hyperplanes in the 
standard way [14-15]; indeed, one can easily derive the usual chessboard estimates [14]. 

In fact A could be built in a similar way out of D-dimensional (hyper) cubes instead 
of squares. Cutting lines become cutting (D — l)-dimensional hyperplanes, and reflection 
symmetry is generalized in an obvious way. Our theorem will then state that the optimal 
flux in each two-dimensional square plaquette of the (hyper) cubic lattice is it in every 
plaquette cut by the hyperplanes. As in the D = 2 case, flux 7r will be optimal in every 
plaquette if we have periodicity in D — 1 directions. 

Returning to the two-dimensional situation (with obvious generalization to D > 2) we 
require only that the U x : s in W° be constant in the horizontal direction (as are the | tail's 
on vertical edges). As for w xy in W d and W s , they are required to be reflection positive 
for reflection in vertical planes (the dashed lines, P, in Fig. 1), as explained below. The 
inclusion of W d and W s is mainly for completeness and nothing essential will be lost by 
setting W d = W s = 0. 

In general, a Hamiltonian, H, can be written (with respect to a cutting hyperplane, 
P) as 

H = H L +H R + H int (5) 

where H L is all the terms involving sites in the left half-cylinder, H R involves right half- 
cylinder sites and H int involves both. The subscripts L, R, and int will also be used for the 
separate pieces, e.g. Kl, Kr, K- mt , etc. Associated with the left Hamiltonian Hl is a right 
Hamiltonian, Q(Hl), which is obtained from Hl by three steps: the unitary transformation 
1Z generated by geometric reflection through the plane P; hole-particle transformation t; 
and complex conjugation *. Q(Hl) = (tIZ^Hl)^ 1 )* . This notion of operator reflection 
can be applied to any operator, Al in the left algebra (i.e., is a polynomial in the 
operators c xa ^c\ (J with x in the left half-cylinder). In particular 0(c la ) = c\ a . A similar 
definition holds for the left Hamiltonian <d(Hn) and, clearly, 0(0(i7^)) = Hl- Note, from 
(4), that S(K L ) = -TI(Kl), Q(W a ) = TZ(W a ), a = 0,d, s. 

Reflection positivity of w xy means that it is symmetric under reflections and, for x 
in the left and y in the right half-cylinder, w xy can be written as a sum (or integral) 
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of functions of the form a(x)a(1Zy)* . In other words, W^ t ,W^ nt are sums (or integral) of 
operators of the form —A L Q(A L ) : where the (— ) in W-^ t comes from r(n x -\ +n x | — l)r _1 = 
1 _ rixi ~ n i|? an d similarly for W? nt . An important example is w xy = w > if x, y are 
nearest neighbors, w xy = otherwise. Such a W s is antiferromagnetic. See [14]. 

Concerning K- int , we note that it is generally not invariant under the three operations. 
However, with / and r denoting a generic left, right image pair cut by P, we are at liberty 
to choose 0(7, r) = 0, i.e., t\ r = |£/ r | > 0, and we do so. (Note: to simplify the notation the 
symbols, j and a, will not be indicated.) This is so because a simple gauge transformation 
c r — > exp[— iO(l, r)]c r makes t\ r > without changing any fluxes. No circuits are involved. 
This choice of phase for ti r is only a convention, for it does not change any physics, but it 
is important for (6) and (7) below. 

With the foregoing convention for H int , the Hamiltonian is said to be reflection sym- 
metric if Q{Hi) = H R . The flux it theorem will be a corollary of the following lemma. 

LEMMA (Reflection positivity). With H as given in (5) with respect to some 
hyperplane, P , assume that K- int satisfies the above positivity convention. Assume also that 
W^ t and Wf nt are reflection positive. Then, for each (3 > and with H int fixed, 

Z(H L ,H R ) 2 < Z(H L ,0(H L ))Z(0(H R ),H R ), (6) 

where Z(Hl, H R ) = Trexp[— j3H]. Moreover, if H R = Q(Hl) and if A^ (resp. A R ) is 
any even operator in the left (resp. right) algebra (e.g., Al is a sum of monomials in cf a 
of even degree) then 

\TrA L A R e-P H \ 2 <TrA L G(A L )e-P H TrG(A R )A R e-P H . (7) 

Proof: Use the Lie- Trotter formula to approximate e~@ H as a product of M > 1 
factors V = V int V L V R , i.e., e~^ H = lim M ^oo V M , where V int = (1 - f3H int /M), V L = 
exp[— /3Hl/M], V r = exp[—(3H R /M]. Notice that Vl contains only even polynomials 
in the c^'s, and so Vl commutes with every right operator (including odd operators). 
Likewise, V R commutes with all left operators. 

If the M factors of Vi nt are multiplied out we obtain for V M a sum of terms, each 
having the form X = a\ Vl V R a,2 Vl V R as Vl V r ••• o,mVl,V r and each ai has one of three 
forms: (i) Al@(Al), with Al an even operator or (ii) c\c r or (iii) — c x c\. Our strategy 
is to move all the left operators to the left without changing the order either of the left 
operators among themselves or the right operators. The operators Al commute with all 
the right operators and cause no difficulty. The difficult point is that the cf operators 
have to move through the cf operators to their left, and each such move gives rise to a 
— 1 factor. 
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I claim that either TrX = or else the number of —1 factors is even. To see this 
note that by particle conservation (and the particle conserving nature of Vl and Vr) the 
number of c\c r factors must equal the number of c x c\ factors if TrX 7^ 0. Call this 
common number J. The number of —1 factors is independent of the order of these 2J 
factors and their order relative to the Al@(Al) factors. The first cf must move through 
zero c^f's. The second cf moves through one cf, etc.. Thus, the number of —1 factors 
is0 + l + 2 + -- - + (2 J — 1) = J(2 J — 1). On the other hand, each c z c\ term carries a 
— 1 factor and there are J of these. Altogether there are J + J(2 J - 1) = 2 J 2 = 0(mod2) 
factors of —1, as claimed. 

In brief, X can be brought into the form X = X^Xr with Xl and Xr even opera- 
tors. Since Trl = 4l A l, we have 4l A lTrA = TrA L TrA fl . Moreover, (TrX L )* = Tr6(A L ) 
and thus (TrA^ 2 = TrX l@(Xl)- Now, denoting the various A's by A a , we have 

\TrV M \ 2 = |£ a TrX Q | 2 = 4_2|A 'I J2 a TTX L TtX R\ 2 < 4 " 2 ' A| E a l Tr ^L| 2 Eal TrX fl| 2 = 

ETrX£6(X£)E a TrX-e(X-) - Z(H L ,0(H L ))Z(H R ,0(H R )). (7) is obtained in the 

Of 

same way. QED 

THEOREM (Flux tv is optimal). Assume the are reflection invariant w.r.t. 
P. Assume also that 0(W") = and W^ t is reflection positive, a = 0, d, s. Then Z is 
maximized by putting flux it in each square face of A that intersects P. 

Proof: We make the gauge transformation above so that K int has tw = \ti r \. From 
(6), we have that when H L ,H R is optimal, so is H L ,0(H L ) and Q(Hr), Hr. But the 
statement K R = Q(K L ) implies the flux it condition by (4). QED 

Remarks: (i). It is interesting to note that (6) can also be used to show that when 
the fluxes are fixed at it and one varies over the \t xy \, the lowest energy is attained in a 
reflection symmetric configuration of \t xy \. A one-dimensional version of the lemma was 
proved and used in [16] to study the Peierls instability for the Hubbard model on a ring. 

(ii) . The lemma and theorem say that flux 7r for and is optimal. If we fix 
<f)*(x,y), we are then free to choose <J)*(x,y) = (f)'(x,y), since any other choice with flux 
7T differs from <f>^(x,y) by a trivial gauge transformation, c x y — > e l ^ x 'c x i. Thus, when 
\txy \ ~ \txy\i ^ ne minimizer can have = T* , thereby preserving SU(2) invariance. 

(iii) . To discuss the ground state we let (3 — * 00. Do we get N = |A| or does the g.s. 
belong to N = |A| + m and N = |A| — m with m > 0? In the Falicov-Kimball model, 
generally, the g.s. has N = 2\A\ and N = 2\B\, as the only choices [17] (but note that 
2\A\ = 2\B>\ = |A| in our case). In the following cases I can also prove that at least one 
g.s. has N = |A|. 

First, assume reflection symmetry and positivity w.r.t. all hyperplanes parallel to P, 
so that we are now looking at a AT with flux it in every plaquette of A. After a trivial gauge 
transformation, this condition can be realized with real T, which we assume henceforth. 
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Next, assume T T = T l , so that SU(2) invariance holds. Third, assume W d = W s = 0, 
i.e., the Hubbard model. As is well known, we can then construct another set of SU{2) 
generators — the pseudospin. (See [18] and, for more details, [19].) By using the spin and 
pseudospin raising operators J2 X c x-[ c x^ ^2x(~^) Xc xt c xi^ (with (-l) x = +1 for x G A, —1 
for x G B) and their adjoints, one can conclude that the absolute ground state belongs to 
A = |A| or A = |A| ± 1. 

Finally, to show that the g.s. has N = |A|, assume either that U x < for all x 
or U x > for all x. We can then use spin-space reflection positivity [17] in Fock space, 
together with the evenness of |A| in our case, to infer N = |A|. (This reflection positivity 
tells us that if a g.s. has numbers iV-j- = A, Aj_ = \i then there are ground states with (A, A) 
and (fx, fx). Thus, if |A| = 2m and N = |A| — 1, so that = m, iVj = m — 1, then there 
is also an (m,m) g.s..) If all U x ^ 0, the g.s. is unique [18]. 

Extensions: The flux phase for the half-filled band has been proved here for a large 
class of Hamiltonians, including the ones common in the physics literature. The proof is 
sufficiently simple that it obviously applies to many other models. 

One generalization is to fermions with n^2 colors, i.e. from SU(2) to SU(n). 

Certain specialized forms of electron-phonon interactions can be included. 

Another generalization is to SU(2) instead of £7(1) gauge fields [15, 20]. Thus, 
^E/Lv is ^placed by \t xy \J2ar c t<r U xy c yr with \t xy \ given, as before, and with 
U xy G SU{2) to be determined. (Even more generally, we can replace t^yC^c ^ +ti y c li c yi 
by E o\t c xa^-x\xy\^xy \^-y\xy c yr where M x - xy G SU(2) and \t xy \ is given; in this case the 
SU(2) matrix associated with xy is U xy = M x . xy M y . xy .) Again, we will find that the en- 
ergy is minimized by flux n in each plaquette, i.e., the product of the four matrices around 
a plaquette satisfies U xy U yz U zw U wx = — 1. 

Thanks are due to I. Affleck, V. Bach, J. Bellissard, E. Carlen, J. Frohlich, M. Loss, 
J.B. Marston, B. Nachtergaele, J. P. Solovej and P. Wiegmann for helpful discussions and 
to the U.S. National Science Foundation, grant PHY90-19433-A03, for partial support. 
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FIGURE CAPTION 
FIG. 1. Typical 2D lattice with horizontal periodic boundary conditions (left boxes 
= right boxes) . Different bond weights illustrate the requirement of horizontal periodicity 
2. Dashed lines (P) are a reflection plane. A generic left-right pair of sites is indicated by 
Z,r. 
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